Landau Quantization in Twisted Bilayer Graphenes: the Dirac Comb 
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We study the Landau quantization of the electronic spectrum for graphene bilayers that are 
rotationally faulted to produce periodic superlattices. Commensurate twisted bilayers exist in two 
families distinguished by their sublattice exchange parity. We show that these two families exhibit 
distinct Landau quantized spectra distinguished both by the interlayer coupling of their zero modes 
and by an amplitude modulation of their spectra at energies above their low energy interlayer 
coherence scales. These modulations can provide a powerful experimental probe of the magnitude 
of a weak coherence splitting in a bilayer and its low energy mass structure. 
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The massless Dirac model that describes the low en- 
ergy electronic physics of single layer graphene is gen- 
erally preempted by the interlayer motion of the charge 
carriers in few layer graphenes (FLG's) This physics 
is well appreciated for Bernal (AB) stacked bilayer 
graphenes where the interlayer coherence scale is ~ 0.2 
cV and the low energy electronic physics 043 is readily 
distinguishable from that of single layer graphene @, @]. 
Surprisingly, one finds only relatively weak (if any) ef- 
fects of coherent interlayer motion in rotationally faulted 
(twisted) FLG's where the crystallographic axes in neigh- 
boring layers are rotated by angles 9 ^ nir/3 [§- 13 1. The- 
ory suggests that the interlayer coherence scale is sup- 
pressed in these structures but below this scale 
one confronts a commensuration problem where the elec- 
tronic physics is controlled both by the rotation angle and 
the atomic registry of neighboring twisted layers 3, li| . 
The low energy electronic properties of these structures 
are very rich but they are not yet well understood 12 , 2(| ■ 



In this Letter we show that the interlayer coupling in 
twisted bilayers produces novel fingerprints in their Lan- 
dau level spectra. We show that small differences in 
the Landau level spectra from two overlapping weakly 
coherence- split bands generate an interference pattern 
detectable as an amplitude modulation of the spectrum: 
the "Dirac comb." The period of the modulation can 
greatly exceed the underlying (presumably small) coher- 
ence scale and its phase contains information about the 
low energy mass structure. This opens the possibility 
of probing even weak interlayer coherence in these sys- 
tems using spectroscopy above their interlayer coherence 
scales. 

Landau quantization of a twisted bilayer is studied us- 
ing a long wavelength theory that includes the effects 



of the lattice commensurability 18]. We study com- 
mensurate rotational faults where we fix overlapping A 
sublattice sites at the origin of the AA stacked stucture 
and rotate one layer through angle 9 with respect to 
the other, generating a new commensurate supcrlatticc. 



Commensurate rotations can be classified according to 
their sublattice parity: sublattice even (SE) structures 
are invariant under exchange of sublattice labels A and 
B while sublattice odd (SO) structures break this sublat- 
tice symmetry. In this notation AA stacking (all sites in 
neighboring layers eclipsed) has 9 = (SE), Bernal AB 
stacking has 9 = n/3 (SO). Sublattice parity can be used 
to classify any commensurate rotated structure and de- 
termines the form of the allowed momentum-conserving 
interlayer terms in the low energy Hamiltonian [Til ] . Each 
layer has Dirac cones in two valleys centered on the K 
and K 1 points described by the (unrotated) long wave- 
length Hamiltonians [2l[ 
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where the 2 x 2 a matrices act on the sublattice (pseu- 
dospin) degrees of freedom. For the SE structures the 
interlayer motion is valley preserving and is described by 
the interlayer operator 



flg? = Ve M exp {icpa z r z 



(2) 



where the 2 x 2 r matrices act on the valley indices. 
For the SO structures the coupling is both interlayer and 
intervalley and is represented by the interlayer matrix 
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(or its valley reversed partner a x Hf^a x ). In Eqns. (2) 
and (3) the coupling strength V, and the pseudospin ro- 
tation angle cj> are determined by the details of the mi- 
croscopic Hamiltonian while the overall phase i9 can be 
removed by a U(l) gauge transformation. 

The low energy dispersion relations for these two fam- 
ilies are shown in the inset of Figure [TJ The spectra of 
SO bilayers have a pair of bands gapped by the interlayer 
coupling and a contact point between two quadratically 
dispersing low energy bands, analogous to the situation 
for the Bernal stacked bilayer except for a ^-dependent 
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FIG. 1: Inset: Dispersion of the coherence-split low energy 
bands for commensurate twisted graphene bilayers that are 
even (blue, solid) and odd (red, dashed) under sublattice ex- 
change. The plot gives the dimensionless energy e — E/V 
as a function of the scaled momentum Q — hvFq/V where V 
is the interlayer coherence scale and the pseudospin rotation 
parameter <f> = 7r/4 for the SE structure. Main figure: Lan- 
dau level spectra Eqns. @ and (|10p for twisted bilayers for 
SE and SO families as a function of the field strength with 
coherence scale V = lOmeV. The red lines are states spanned 
by the zero modes of the two layers. 



q = with equal weights on the two layers and a coher- 
ence splitting produced by their interlayer coupling. Note 
that the SE bilayers are generically gapped systems (their 
band extrema occur on a ring in reciprocal space) due to 
an avoided crossing of its two coherence-split branches. 
The size of this gap is determined by the pseudospin ro- 
tation angle <f> in Eqn. (2) and it vanishes for the special 
case of AA stacking where <j> = by symmetry. 

To study this system in a perpendicular magnetic field 
B we introduce a vector potential in the symmetric 
gauge A = (B x r)/2. Defining z = (x + iy)/lB with 
magnetic length £b = \JhJeB and cyclotron frequency 
uj c = \^2vf/£b the Hamiltonian at the K point of the j- 
th layer (j = 1,2) is transformed to a Landau level (LL) 
basis 



h,K,j = —ioj c 
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where at — (—2d z + z/2) /-y/2 is the Landau level rais- 
ing operator and 8j is the rotation angle of the j-th 
layer. The corresponding Hamiltonian at the K' point 
is hx'.j — CTyfiKjCTy Since the interlayer coupling matri- 
ces in Eqns. (2) and (3) are local in the layer-projected 
coordinates the bilayer LL Hamiltonian has the form 
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for SE bilayers and in the SO case 
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interlayer coherence scale V which can be significantly re- 
duced. The SE structures feature a pair of Dirac nodes at 



Particle-hole symmetry allows the spectrum to be studied 
by squaring Ji, yielding for SE bilayers 
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where evi — Q\ +cj> and 02 = #2 — 4> an d for the SO bilayer 
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along with their valley reversed partners. Eqns. and 
© demonstrate that for either family of structures the 
Hamiltonian is block diagonal in a particular layer- and 
valley-polarized Landau level basis. For SE structures 
the interlayer coupling scatters states from LL n in one 
layer only into states from LL±„ in the same valley of 
the neighboring layer, while in the SO structure LL„ in 
valley K of layer 1 couple only to LL±( n _i) in valley K' 



of layer 2. Within these subspaces the Hamiltonians © 
and (J5J can be diagonalized yielding the spectra 

el B KtU = ±\jnuj 2 c + V 2 + nV^/2nuj 2 c {\ + cos0) (9) 

and 
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with valley (y = ±1) and branch [k — ±1) indices, set- 
ting H = 1. The eigenstates of the SE structures come in 
valley-degenerate pairs. SO structures have degenerate 
pairs of states along with a quartet of orbital states at 
e = (with indices n = 0, k = — 1, and v = ±1). 

In Figure [T] we plot the LL spectra of Eqns. © and 
(fT0| as a function of the held strength. For SE structures, 
the states disperse oc \[B away from the two coherence- 
split Dirac nodes and evolve at high field into two pairs of 
branches that disperse away from a common charge neu- 
trality level. For SO structures the fermions are massive 
in both the low energy (around E = 0) and high energy 
(around E = ±V) branches and the LL dispersion is thus 
linear in B when uj c <C V crossing over to the expected 
the \B dependence in high held. For Bernal stacking 

V > w c at experimentally realizable fields and one can 
integrate out the high energy band to obtain an effective 
two band model for the "weak field" limit of the SO spec- 
tra 0. By contrast in twisted SE and SO structures the 
coherence scale collapses and all the degrees of freedom 
become accessible. 

A striking prediction of Eqns. ((9]) and (fP0|) is that 
the overlapping branches of Landau levels produce an 
amplitude modulation of the density of states. When e ^> 

V for SE structures and e 3> (V, w 2 /V) for SO structures 
we can approximate e ps ±(y/n+n(3V) where for (SE,SO) 
bilayers (3 — (cos(</>/2), 1/2). The density of states is 
enhanced whenever the energy of LL„ in one branch (k — 
— 1) overlaps the energy of LL„_ m in the other (k = 
1). This occurs for m = AfiV^/n/Lo c when n 3> m, and 
produces a beating pattern with period AE = w 2 /4/3V, 
most evident when AE ui c 3> V. The superposition of 
the LL spectra thus produces a "Dirac comb." 

Scanning tunneling spectroscopy (STS) [22] in mod- 
est magnetic fields can access these quantum oscilla- 
tions in the single layer-projected density of states p\ 
and thus quantify the interlayer coherence. Theoret- 
ically, p\ is conveniently studied by integrating out 
the second layer [23| yielding an effective Hamilto- 



nian for the first (exposed) layer Hf (E) = Hi + 
H m t (E — Hz) -1 H- lnt containing a state- and energy- 
dependent self energy that is evaluated separately in 
each invariant subspace of Eqns. (0 and (j6]). The 
density of states is then obtained by a trace over the 
sublattice degrees of freedom and invariant subspaces: 
Pi(E) = E„,„Imtr [E - ij - Hf% >u (E - i 7 )] _1 with 
level broadening 7. Importantly, in both SE and SO 
structures the zero mode of the surface-projected prob- 
lem occupies a one dimensional subspace with 

|V|7£; Hl s n i ^ 8 V ^E\V\ 2 /{E 2 -Lol).{ll) 

Eqn. (fTTj) demonstrates that the zero modes of each layer 
are always coherence-split by 2V in SE bilayers but re- 
main exactly decoupled for all SO bilayers. 

These features can be seen clearly in the densities of 
states in Figures (2) and (3). The coherence splitting 
of the zero modes for the SE structure can be contrasted 
with the degenerate layer-polarized zero energy states for 
the SO structures. At higher energy (Fig. (3)) quantum 
oscillations are clearly seen for both families. Although 
the B = dispersions are nearly the same for these two 
structures at energies E S> V (Fig. 1) the amplitude 
modulations near E ~ oj 2 //3V are phase shifted, reflect- 
ing their distinct low energy mass structure. 

The modulations reflect the atomic registry in a com- 
mensurate bilayer and are most pronounced for short 
period superlattices. They are observable for B > 
V 2 /2ehvp, i.e. at a field scale of only 0.1 T for V ~ 
lOmeV. For small fault angles the coherence scale V 
collapses, the magnetic field scale is correspondingly re- 
duced and measurements will ultimately be limited by 
the finite quasiparticle lifetime. Interestingly, the mod- 
ulation period AE = h 2 cu 2 /4(3V = HeB/M* describes 
the cyclotron frequency of a massive particle with M* ~ 
3.5 x 10~ 3 77i c when V = lOmeV. This can be understood 
semiclassically by noting that the quantization conditions 
for orbits in the two weakly coherence-split Dirac bands 
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FIG. 2: Layer projected density of states for SE (blue) and 
SO(dashed red) twisted bilayers as a function of E /uj c (h = 1). 
Data are plotted for V /uj c = 0.1, 7 = .03 x uj c and (f> = 7i"/4. 
The low field states of the SO bilayer are weakly coherence 
split for the SE bilayer. 




gularity where the topology of the bands changes; in a 
magnetic field this is identified by an onset of enhanced 
coherence-splittings of nearly layer-degenerate LL states. 
For a rotation angle ~ 1°, AE ^150 meV so for the n-th 
Landau level this occurs for nB ~ 30 T. Likewise the 
physics of the Dirac comb is distinct from the Hofstadter 
physics that arises when the magnetic length becomes 
commensurate with a superlattice translation, requiring 
B ~ 4Tfor 6 ~ 1° [H. 

In summary we have studied the Dirac comb in the 
weak field Landau quantized spectra for twisted graphene 
bilayers: an interference phenomenon yielding an ampli- 
tude modulation of the LL spectra. This is important fol- 
low energy magnetotransport and can be used to iden- 
tify weak interlayer coherence and its low energy mass 
structure. 

This work was supported by the Department of Energy, 
Office of Basic Energy Sciences under contract DE-FG02- 
ER45118 (EJM) and by the NSF DMR-0820382 (MK). 




FIG. 3: Layer projected density of states as a function of E/u) c 
over a wide energy range showing the amplitude modulation 
in the Dirac comb. Data are plotted for V /ui c = 0.3 and 
<f> = 2tv/3. Top panel(blue, vertically offset) is for an SE 
bilayer and the lower panel (red) is for an SO bilayer. 



acquire the difference of one quantum from one maximum 
of the amplitude modulation to the next. An effective 
quantization condition within the phase space annulus 
bounded by the two coherence-split bands thus follows. 

These low field quantum oscillations are physically dis- 
tinct from the Landau level structure studied for twisted 



graphene bilayers in high fields |24| . Low angle faults off- 
set the Dirac nodes at the zone corners K in neighboring 
layers introducing an energy scale AE = HupK sin(0/2) 
at which layer-decoupled Dirac cones intersect. Inter- 
layer mixing of these states creates a saddle point sin- 
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